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Abstract 

A kinetic theory for quantum Langmuir waves interacting nonlinearly with quantum ion-acoustic 
waves is derived. The formulation allows for a statistical analysis of the quantum correction to the 
Zakharov system. The influence of a background random phase on the modulational instability 
is given. In the coherent case, the effect of the quantum correction is to reduce the growth rate. 
Moreover, in the classical limit, a bifurcation develops in the dispersion curves due to the presence 
of partial coherence. However, the combined effect of partial coherence and a quantum correction 
may give rise to an increased modulational instability growth rate, as compared to the classical 
case. The results may be of significance in dense astrophysical plasmas and laboratory laser-plasma 
systems. 

PACS numbers: 52.35.-g, 03.65.-w, 05.30.-d, 05.60.Gg 
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I. INTRODUCTION 



The effects from the quantum domain have intriguing consequences for the way we view 
the world and the way we interpret physical models, and these effects can be seen in single 
particle phenomena. Plasma physics, on the other hand, deals with the collective interaction 
of charged particles, and processes related to such. However, there are physical parameter 
domains where quantum mechanics and plasma physics need to be taken into account simul- 
taneously, e.g. dense and/or hot astrophysical plasmas. Indeed, the creation of pair plasmas 
surrounding neutron stars can be viewed as a collective quantum plasma effect. Thus, there 
is interest and application for a model taking both collective charged particle effects and 
quantum phenomena into account (see, e.g. Refs.QandS). 

One of the most prominent models in plasma physics is described by the Zakharov 
equations^, in which high frequency Langmuir waves are coupled nonlinearly to low fre- 
quency ion-acoustic waves. The statistical properties of this system has been analyzed in 



Ref. U, where a Landau-like damping was found. Recently, a generalization of the Zakharov 
system was derived, taking quantum effects into account^. The effect of the quantum cor- 
rection was to introduce higher order dispersion into the system of equations, thus altering 
the behavior of wave evolution. It was arg 

important in astrophysical plasmas, as the plasma densities may become significant. 

In this paper, we will introduce a kinetic description of the quantum Zakharov equation, 
by applying the Wigner transform to the Langmuir propagation equation. The resulting 
system of equations may be useful for understanding the properties of partially coherent 
Langmuir waves interacting with quantum ion-acoustic waves. We derive the general dis- 
persion relation, and analyze the stability of the system. A comparison between monoener- 
getic Langmuir waves and random-phase Langmuir waves is given, and it is found that the 
interplay between quantum corrections and spectral broadening may alter the instability 
properties in novel ways. In particular, it is found that the growth rate for a short wave- 
length partially coherent quantum Langmuir wave is larger than the corresponding growth 
rate for the classical Langmuir wave. Thus, the interplay between incoherence and quan- 
tum effects gives rise to modified modulational instability growth rates, a result that may 
be relevant to astrophysical and intense laboratory laser-plasmas, for which the quantum 
parameter may take on significant values. 
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II. BASIC EQUATIONS 



The dynamics of the nonlinearly coupled quantum Langmuir and ion-acoustic waves is 
given by the Zakharov-like equations 5 

id t E(t, x) + d 2 x E{t, x) - H 2 d*E{t, x) = n(t, x)E(t, x), (la) 

and 

(d 2 - %)n(t,z) + H 2 din(t,x) = d 2 x \E(t,x)\ 2 } (lb) 

where H = hu p i/kf,T e is the quantum parameter due to a quantum pressure emanating 
from the underlying hydrodynamic model 5 ". Here % is Planck's constant divided by 2tt, 
u P i = (noe 2 / rriieo) 1 ^ 2 is the ion plasma frequency, fee is Boltzmann's constant, T e is the 
electron temperature, no is the constant background density, and is the ion rest mass. 
The electric field E has been normalized according to E — > {e^rrii /IQrrienokBTe) 1 / 2 E , while 
the density n is normalized by n — > (m i /4m e no)n, where m e is the electron mass. The 
coordinates have been rescaled using t — > (2m e /mi)uj pe t and x 2(m e /mj) 1//2 x/A e , where 
ujpe = (n,Qe 2 /m e eo) 1//2 is the electron plasma frequency and A e is the electron Debye length. As 
H — > 0, we regain the classical Zakharov equations from ((TJ. However, in some astrophysical 
plasmas, the quantum parameter H may approach unity, since in such environments, high 
densities are not uncommon (see, e.g. Ref. 0). We see that the effect of the quantum 
parameter is to introduce higher order dispersion. 



III. QUANTUM KINETICS 

The Fourier transform of the two-point correlation function, as given for the electric field 

by 

p(t,x,p) = ± J d£e*(Er(t,x + t/2)E{t,x-t/2)) (2) 

was introduced by Wigner«£ in quantum statistical mechanics. Here the angular brackets 
denotes the ensemble average, and the asterisk denotes the complex conjugation operation. 
The Wigner function p is a generalized distribution function, which satisfies 

{\E(t,x)\ 2 )= [ dpp(t,x,p). (3) 
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Applying the transformation (J2J) to Eq. ()la|) gives the kinetic equation 8 

dtp(t, x, p) + (2pd x + AH 2 p 3 d x - H 2 pdl) pit, x, p) - 2n(t, x) sin (j dj^j pit, x, p) = 0, (4) 

which is coupled to the ion-acoustic equation (fTb"j) via Eq. ©. Here the sin-operator is 
defined by its Taylor expansion, and arrows denote direction of operation. Keeping the 
lowest order derivative in this Taylor expansion, corresponding to the long wavelength limit, 
gives a modified Vlasov equation 

d t p(t,x,p) + pD x p(t,x,p) - {d x n{t,x)){d p p{t,x,p)) = 0, (5) 

for the quantum Langmuir wave, driven by the ion-acoustic ponderomotive force. Here 
D x = (2 + AH 2 p 2 — H 2 d 2 )d x . Thus, in the classical limit H —>■ 0, D x —>■ 2d x , and we obtain 
a Vlasov-like equation for the long wavelength Langmuir waves. 



IV. THE MODULATIONAL INSTABILITY 

In order to analyze Eqs. (jlb|) . ©, and (@J), we perform a perturbative expansion. Letting 
p(t,x,p) = Po(p)+Pi exp(ikx— iut), where \pi \ <C po, and n(t, x) = n +rii exp(ikx — iut), we 
linearize with respect to the perturbation variables. We then obtain the dispersion relation 

The dispel .elation Q g ™^» the results m Rets.OandB and is valid for partially 
coherent quantum Langmuir waves interacting nonlinearly with quantum ion-acoustic waves. 



A. Monoenergetic Langmuir waves 



In the case of a monoenergetic Langmuir wave, we have po(p) = Iq5(p),, where Iq = \E \ 2 
is the background intensity. Then the dispersion relation (JBJ) becomes 

[oo 2 - (1 + H 2 k 2 )k 2 ] [u 2 - (1 + H 2 k 2 ) 2 k 4 ] = 2/ (l + H 2 k 2 )k\ (7) 



such that 



2 _ 1 



H 2 k 2 + H A k A ± Hk 2 \/H 2 + 8/ - 2H 4 k 2 + H e k* 



(8) 
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where H = 1 + H k . Letting to = ry, the instability growth rate is given by- 



7" 1 



v/2 



Hk 2 ^H 2 + SI - 2H 4 k 2 + H Q k 4 - H 2 k 2 - H^k 4 ] ^ . (9) 



Starting from H = 0, successively higher values of H tend to suppress the instability, giving 
lower growth rates with a cut-off at a lower wavenumber, see Fig. ^ 

B. Partial coherence 

The coherent monoenergetic background distribution gives important information on 
wave instabilities. However, in many applications the background field is not fully coherent, 
but rather displays partial decoherence due to, e.g. noise. The noise, either classical or quan- 
tum, may stem from different sources, such as thermal effects, weak turbulence, or quantum 
fluctuations. Such sources of noise may lead to a background field E with a random phase 
ip(x) such that 

( e -i[<p{x+t/2)-<p(x-£/2)]j _ e -p w \t\^ (10) 

with the corresponding distribution function po is given by the Lorentzian&iJi 

P0{P) = 2 1 2 > 11 

where pw is the width of the distribution. The integrand of © has three poles, one real 
and two complex, where the real pole is given by 

1 A 1/3 (uj, k) 

P °-~A^(u,k) + 3b(k) ■ (12) 

Here A(u, k) = 3[9ab 2 + v / 3(46 3 + 27a 2 6 4 ) 1/2 ]/2, a(u, k) = u/k{2 + H 2 k 2 ), and b{k) = 
AH 2 / (2 + H 2 k 2 ). As the quantum parameter H approaches zero, the complex poles ap- 
proaches complex infinity. Thus, in the integration of Eq. (jUJ) we will neglect these poles, 
only taking the real pole po into account, since the modes corresponding to the complex 
poles are quickly damped. Thus, we have 



-a/ + (1 + H 2 k 2 )k 2 = 2I k 



g{k)-h{k) 



[u + 2ip w h(k)k] 2 - k A [g(k) - h(k)f 
ipwPo 



k\2 + H 2 k 2 )[(p + k/2) 2 +p w ][(p - k/2) 2 + Pw ] 



(13) 
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where we have defined the real and positive functions g(k) = H 2 (k 2 + 8p 2 v ) and h(k) = 
1 + 2H 2 (k 2 +p 2 v )- The dispersion relation (|TB"|) describes the effects of partial coherence for 
the quantum Zakharov system (JTJ). The damping character due to the finite width of the 
background distribution can clearly be seen, as well as the Landau damping due to the real 
pole. We note that as pw — ► 0, we regain the monoenergetic dispersion relation (J7J). 

C. The classical limit 

If H — > 0, we obtain the classical limit of the dispersion relation (|13|) . when we use a 
kinetic photon description for the Langmuir waves. The effects of statistical broadening 
on the Zakharov system was also analyzed in Ref. |4J We note that the two complex poles 
approaches infinity, and only the real pole remains with the value po = uu/2k, as it should. 
The dispersion relation then reads^ 

(to 2 -k 2 ) [u 4 + k 8 + 8p 2 w k 6 + 8p 2 w k 2 u 2 + 2k\8p 4 w - u 2 )] = 2I k\u 2 -k 4 -8ip w ku-Ap 2 w k 2 ) ) 

(14) 

and although the quantum effects have been neglected, the behavior of the function u(k;pw) 
is still rather complicated. The dispersion relation (0) with H = was analyzed analytically 
in the long wavelength limit, i.e. u/k ^> 1, in Ref. y, and the growth rate was found for a 
Gaussian background spectrum. Here we will solve the equation ()14j) for all wave lengths. 
Letting uo = Re cu+27, we may solve Eq. (|13|) numerically for the growth rate 7, using different 
values of the widths pw- In the Figs. El and El we have plotted the solutions for H = and a 
number of values of pw The results show on a more complicated dispersion structure than 
in the coherent case. The asymptotic behavior of the growth rate for short wavelengths has 
been depicted in Fig. HJ using a number of different values on the decoherence width pw- 
For large k, the growth rate has a linear dependence on the wavenumber, the slope being 
determined by the values of the width p w . 

D. Quantum effects on the instability growth rate 

When H is nonzero, the combined effects of quantum correction and decoherence make 
themselves apparent in the dispersion relation (|13j) through new and novel wave modes. In 
Fig. (0) we display the growth rate 7 as a function of the wave number k. As compared to 
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the classical case, the combined effect of partial coherence and quantum effects, i.e. finite 
Pw and H respectively, is to make the modulational instability growth rate smaller for long 
wavelength perturbations. However, the interesting effect is for short wavelengths, where 
the modes introduced due the finite spectral width is amplified by the quantum corrections. 
Thus, we may expect much stronger growth rates for short wavelength perturbations, making 
these dominant in quantum plasmas. In Fig. |3] the dispersion curves for a value H = 0.25 
of the quantum parameter has been plotted. The strong growth rate for short wavelengths 
can clearly be seen. 

V. CONCLUSIONS 

The effects of partial coherence in quantum plasmas, such as in the form of a ran- 
dom phase, is of interest in certain plasmas, such as astrophysical plasmas*^, and the 
next generation laser plasma system a 11 ' 12 . Moreover, in such system, the density may 
even reach values of^ 10 23 — 10 31 m~ 3 for temperatures of the order 10 5 — 10 7 K, giving 
H = huo pi /kBT e ~ 10~ 7 — 1. Thus, the quantum parameter H may attain appreciable val- 
ues, such that the higher order dispersive terms in Eqs. (JIJ become important. Even in 
the cases of small a quantum correction, this effect combined with Langmuir wave decoher- 
ence will yield a strongly growing mode for short wavelengths, and could lead to significant 
changes in extreme astrophysical and laboratory plasmas. Thus, the combination of inco- 
herence and quantum effects may yield rich and interesting dynamics of Langmuir wave 
propagation in such plasmas. However, a detailed analysis of possible applications is left for 
future research. 

Here we have analyzed the statistical properties of the quantum Zakharov system, giving 
the dynamics of high frequency Langmuir waves in terms of a kinetic equation. This enabled 
the investigation into the effects of partial coherence of the quantum Langmuir wave, in 
particular the implications due to a random phase, and it was found that such a system 
exhibits an interesting dispersion structure. In particular, the combined effect of decoherence 
and quantum corrections gives rise to new dispersion curves as well as increased modulational 
instability growth rates, as compared to the case of a classical coherent and partial coherent 
Langmuir wave. 
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Fig. 1: The growth rate 7 as given by Eq. Q plotted as a function of the wavenumber k. 
Seen from the right, the values H = 0, 0.5 and 1 have been used together with the intensity 
I = 0.5. The effect of the quantum parameter is thus to suppress the instability. 

Fig. 2: The classical limit of the dispersion relation (fT3*j) for different values of the 
background distribution width p w . We see that the upper curve in panel (a) closely 
resembles the uppermost curve in Fig. ^ but is slightly damped. The effects of the spectral 
width can clearly be seen through the damping of the mode present in the coherent case, 
as well as the presence of a completely new mode. In the different panels, we have the 
following widths: (a) p w = 0.025, (b) p w = 0.05, (c) p w = 0.0667, (d) p w = 0.0909, and 
(e) pw = 0.267. The intensity in all the panels is 1$ = 0.50. 

Fig. 3: The classical limit of the dispersion relation (fT3*j) for a select set of values of the 
background distribution width pw, taken to larger wavenumbers. In the different panels, 
we have the following widths: (a) pw = 0.0714, (b) p w = 0.111, and (c) pw = 0.333. The 
intensity in all the panels is I = 0.50. 

Fig. 4: The inverse of the classical growth rate, i.e. H = 0, plotted as a function of the 
wavelength A = 2ir/k, giving the asymptotic behavior of 7 for short wavelengths. Starting 
from the top of the panel, we have used the values pw = 0.025, pw = 0.05, and pw = 0.091 
on the respective curve. We note the generic linear behavior for short wavelengths. 

Fig. 5 (Color online): The combined effects of a quantum correction and partial coherence 
on the growth rate (black curve), obtained from Eq. (fT3*|) using H = 0.25, as compared to 
the classical case for the same spectral width (red curve). In the different panels, we have 
the following widths: (a) pw = 0.05, (b) pw = 0.111, and (c) pw = 0.2. We note the mode 
due to the quantum corrections combined with the partial coherence gives rise to a larger 
growth rate for short wavelengths, while the long wavelength modes are damped by the 
quantum corrections (see Fig. [TJ). 
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